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ABSTRACT 

Context. What are the mechanisms by which massive stars form? What are the initial conditions for these processes? It is commonly 
assumed that cold and dense Infrared Dark Clouds (IRDCs) likely represent the birth sites massive stars. Therefore, this class of 
objects gets increasing attention, and their analysis offers the opportunity to tackle the above mentioned questions. 
Aims. To enlarge the sample of well-characterised IRDCs in the southern hemisphere, where ALMA will play a major role in the near 
future, we have set up a program to study the gas and dust of southern infrared dark clouds. The present paper aims at characterizing 
the continuuum properties of this sample of IRDCs. 

Methods. We cross-correlated 1.2 mm continuum data from SIMBA@SEST with Spitzer/GLIMPSE images to establish the con- 
nection between emission sources at millimeter wavelengths and the IRDCs we see at 8 /im in absorption against the bright PAH 
background. Analysing the dust emission and extinction leads to a determination of masses and column densities, which are impor- 
tant quantities in characterizing the initial conditions of massive star formation. We also evaluated the limitations of the emission and 
extinction methods. 

Results. The morphology of the 1.2 mm continuum emission is in all cases in close agreement with the mid-infrared extinction. The 
total masses of the IRDCs were found to range from 150 to 1150 M Q (emission data) and from 300 to 1750 M Q (extinction data). 
We derived peak column densities between 0.9 and 4.6 xlO 22 cm 2 (emission data) and 2.1 and 5.4 xlO 22 cnr 2 (extinction data). 
We demonstrate that the extinction method fails for very high extinction values (and column densities) beyond Ay values of roughly 
75 mag according to the Weingartner & Draine (2001) extinction relation Ry = 5.5 model B (around 200 mag when following the 
common Mathis (1990) extinction calibration). By taking the spatial resolution effects into account and restoring the column densities 
derived from the dust emission back to a linear resolution of 0.01 pc, peak column densities of 3 - 19 x 10 23 cnr 2 are obtained, much 
higher than typical values for low-mass cores. 

Conclusions. The derived column densities, taking into account the spatial resolution effects, are beyond the column density threshold 
of 3.0 xlO 23 cnr 2 required by theoretical considerations for massive star formation. We conclude that the values for column densi- 
ties derived for the selected IRDC sample make these objects excellent candidates for objects in the earliest stages of massive star 
formation. 

Key words. ISM: dust, extinction, ISM: clouds, Infrared: ISM, Radio continuum: ISM, Stars: Formation 
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Jjj , 1. Introduction first identified wit h the Infrared Space Observatory (ISO; 

. . . iPerault et al.l [1 996) and Midcourse Space Experiment (MSX; 

One of the key questions in stellar astrophysics is to understand E ga n et alJU998 !) as dark extended features against the bright 

the formation and earliest evolution of high-mass stars. These Galactic PA H backgrou nd at mid-IR (MIR ) wavelengths. The 

objects play a major role in shaping the interstellar medium due first studies dEgan et al.Nl 998k ICarev et aljfl998l I2000 ) showed 

to their strong UV radiation fields and stellar winds, and they that IRDCs ^ dense (>10 5 cm -3 )i cold (<25 K ) and can at- 

ennch their environment with heavy elements when exploding as tain high column densities ( ;> 1Q 23 cm -2 y An these properties 

supernovae. Despite their importance, the mechanism by which make IRDCs excellent candidates for hosting very early stages 

such massive stars form and, especially, t he initial conditions Q f mass j ve star formation 
of their birthplaces are p oorly understood (Beuther et alj|2007t 

IZinnecker & Y orke 2007]). During the last years, additional studies of Infrared Dark 

One pathway to tackle especially the latter problem is to Clouds at millimeter a nd sub millimeter wavelengths were per- 

analyse so-called Infrared Dark Clouds (IRDCs). They were formed. |Simon et al.| (|2006a|) presented a catalog of almost 

1 1 000 IRDCs in the first and fourth quadrants of the Galactic 

Send offprint requests to: T. Vasyunina P lane - Usin g 13c ° J=1 -° molecular line emission, the kine- 

* Based on observations made with the ESO 37-channel bolometer matic distances to 31 3 clouds from this catalog were established 

array SIMBA at the SEST telescope on La Silla, under programme ID dSimon et alj|2006b|). This allowed to estimate sizes, masses and 

71.C-0633. the Galactic distribution for this large sample. The study showed 

** Member of the International Max Planck Research School (IMPRS) that IRDCs have sizes of ~ 5 pc and LTE masses of ~ 5 X 

Heidelberg 10 3 M comparable to cluster-forming molecular clumps. The 
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galactic distribution of IRDCs follows the general distribution 
of molecular gas. A concentration of the clouds is associated 
with the so-called 5 kpc molecular ring, the Galaxy's most mas- 
sive and active star-forming complex. Ammonia observations of 
some well-known IRDCs with the Effelsberg 100 m telescope 
(Pillaj_etaL|l2006a) allowed to estimate additional chemical and 
physical properties, such as average gas temperature (between 
10 and 20 K), velocity fields (significant velocity gradient be- 
tween the cores, linewidths of 0.9-1.5 km/s) and the chemical 
state. According to this study, NH3 in IRDCs is overabundant by 
a factor of 5 - 10 relative to Taurus or Perseus local dark clouds, 
while H2CO is underabundant by a factor of ~ 50. Hence, the 
chemistry governing IRDCs might be complex and could be dif- 
ferent from other parts of the molecular ISM. 

Although significant progress has been made observation- 
ally, the number of IRDCs with well characterised properties is 
still small to date, especially regarding the southern hemisphere. 
To enlarge the sample of well-investigated IRDCs, we selected 
12 clouds in the southern hemisphere and started a program to 
measure the gas and dust properties of these objects. 

In Sect. 2 we describe our source selection and 1 .2 mm con- 
tinuum observations with the SIMBA/SEST telescope. In Sect. 
3, we discuss the data reduction and the details concerning the 
calculation of dust masses and column densities. Also we present 
here the comparison of the MIR and millimeter techniques. In 
Sect. 4, we compare our results with previous results for highl- 
and low-mass star-forming region and with the theoretical mod- 
els. 



2. Observations 

The IRDCs for our study were selected in the "pre-Spitzer" era, 
by visual examination of the MIR images delivered by the MSX 
satellite. The MSX A band (6.8 - 10.8 //m) was the most sensi- 
tive one among the MSX bands and presented the highest level 
of diffuse background emission (due to PAH emission at 7.7 and 
8.7 //m), which also leads to the highest contrast between bright 
background and dark IRDCs. We selected a sample of southern 
IRDCs from the A band images looking for high contrast and 
sizes sufficient to fill the main beam of the SEST telescope at 
1.2 mm. 

The 1.2 mm continuum observations were carried out w ith the 
37-channel bolometer array SIMBA dNvman et alj|200ll) at the 
SEST on La Silla, Chile between July 16-18, 2003. SIMBA is a 
hexagonal array in which the HPBW of a single element is about 
24" and the separation between elements on the sky is 44". The 
observations were mad e using a fast mappin g technique without 
a wobbling secondary (Weferling et al. 2002). 

Maps of Uranus were taken to check the flux calibration of 
the resulting data. To correct for the atmospheric opacity, sky- 
dips were performed every 2-3 hours. Despite the occurrence 
of some thin clouds, the observing conditions were good which 
is reflected in zenith opacity values of 0.16-0.18. The point- 
ing was checked roughly every two hours and proved to be bet- 
ter than 6". The combination of typically three maps with sizes 
of 560" x 900" resulted in a residual noise of about 22-28 
mJy/beam (rms) in the center of the mapped region. 



3. Data reduction and analysis 

In this paper we are using both 8 fim IRAC data from the Spitzer 
Galactic Le gacy Infrared Mid-P lane Survey Extraordinaire 
(GLIMPSE, Benj amin etall 120031) and our 1.2 mm data from 



the SIMBA bolometer at the SEST telescope to investigate the 
physical properties of the extinction and emission material. 

In both cases, for deriving the masses of the IRDCs, we need 
a handle on the distances to these clouds. For determining the 
(kinematic) distances to our IRDCs, we use the vlsr veloci- 
ties derived from recent molecular line observation^] with the 
Australian MOPRA telescope, which we will present in a forth- 
coming paper. The velocities have been transfered to kinematic 
distances by adopting th e recently improved parameters for the 
Galactic rotation curve dLevine et alj [2008) for the fourth and 
first Galactic quadrant. Always the near kinematic distance has 
been assumed. The corresponding distances are reported in Table 
[3]and have been used for the mass estimations. Note that such ro- 
tation curves give just average properties. The actual distribution 
of material might be more structured, especially in the fourth 
quadrant, which is indi c ated in the HI absorption measurements 
shown in iLevine et alj (2008). Furthermore, for objects in the 
Galactic longitude interval [305°, 310.5°], several velocity sys- 
tems can occur due to the projection of at least two Galactic 
arms. The measured velocities, however, place all our IRDCs in 
that longitude r ange within the Ce ntaurus arm (3.5 - 5.5 kpc), in 
agreement with lSaito et alj d2001l) . 

3.1. Millimeter data 

The 1 .2 mm data for the IRDC regions from SIMBA at the SEST 
telescope were reduced using the MOPSI package (developed 
by R. Zylka, IRAM). All maps were reduced by applying the at- 
mospheric opacity corrections, fitting and subtracting a baseline, 
and removing the correlated sky noise. Thereby, we followed a 
three-stage approach as suggested in the SIMBA manual. After 
a first iteration using all data for the sky noise removal, the map 
regions showing source emission are neglected for sky noise re- 
moval in the second iteration. From this second interim map a 
source model is derived which is being included in the third iter- 
ation. The resulting maps were flux-calibrated using the con- 
version factor obtained from observations of Uranus. For our 
July 2003 observations, this factor was around 60 mJy/beam per 
count. 

For estimating cloud masses and column densities we used 
the following expressions: 



M tn , = 



Pint D 2 R 
B V (T) x v ' 



i peak 



R 



Q B V (T) x v mn. 



(1) 



(2) 



^peak is the measured source peak flux density, while F[ nt 
denotes the integrated flux density of the whole source. D. is the 
beam solid angle in steradians, mn, is the mass of one hydrogen 
molecule, D is the distance to the IRDC, R is the gas-to-dust 
ratio, x v is the dust opacity per gram of dust, and B V (T) is the 
Planck function at the dust temperature T. We adopt a gas-to- 
dust mass ratio of 100, and x v equal to 1.0 cm 2 g~', a valu e 
appropriate for cold dense cores (Ossenkop f & Henning| [T994). 

At the present stage, where measured temperatures are not 
available for our sources, we assume a temperature of 20 K, 
which is a reasonabl e choice considering recent investig ations 
toward other IRDCs dCarev et all l998; Pilla i et al.ll2006bi) . 
The derived mass depends on the (assumed) temperature, on the 



The HCO + (1-0) line velocities have been employed for this pur- 



pose. 
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distance to the cloud, and on the grain model. Masses will be un- 
derestimated if the temperature is lower than the assumed value 
of 20 K. For example, for 15 K the masses will be higher by a 
factor of 1.5. In case of a higher temperature in the cloud, e.g., 
for 30 K, our results have to be divided by a factor of 1.7. We 
further note the quadratic dependence of the derived masses on 
the distance of the clouds. Hence, the masses will be by a factor 
of 1.2 - 1.8 higher if the distance is 500 pc more, and lower by 
the same factor if it is 500 pc less than indicated by the average 
Galactic rotation curve (see above). The derived masses are in- 
verse proportional to the assumed value of the opacity x v , which 
has an uncertainty of at least a factor 2. The column density has 
no direct dependence on the distance to the cloud, but the tem- 
perature dependence is the same as for the masses. 



Table 1. Comparison of the 1.2 mm emission and 8 /im absorp- 
tion techniques for deriving masses and column densities. 



Properties 



1.2 mm 



Distance dependence" 
Resolution 

Sensitive to lower column 
density filigree structure 

Temperature dependence 

Background and foreground 
estimation 

Sensitive to column 



+ 

24" 



+ 

not 
necessary 

+ 



8 /im 



+ 
3" 
+ 



necessary 



3.2. GLIMPSE 8 /jm data 



Only for masses 



The original selection of the IRDCs has been done still based 
on MSX images. In the meantime, the Spitzer satellite has suc- 
ceeded MSX and provides much higher spatial resolution and 
sensitivity. GLIMPSE images for our regions with a pixel size 
of 076 were retrieved from the NASA/IPAC Infrared Science 
Archive (IRSA) and re-mosaicked to cover the final field of our 
1 .2 mm maps of typically 15' x 15' in order to get a picture of the 
IRDCs in relation to their closer and further vicinity (see Figs. 
HH2). After bad pixel removal, a PSF phot ometry has been per - 
formed using the STARFINDER program dDiolaiti et alJ feOOO). 
This allows to remove compact foreground objects, and thus 
maps of extended emission and absorption structures and finally 
column density maps could be extracted in a subsequent step. 

Dust masses for the IRDCs were computed assuming that 
they are in the foreground and shadow emission from behind. 
The optical depth t is ideally the logarithm of the ratio of two 
fluxes: (a) the flux from the emission background 7b ac k directly 
behind the IRDC, and (b) the actually measured remnant flux 
/irdc from the location of the IRDC. Furthermore, superimposed 
on the IRDC is an emission contribution from foreground mate- 
rial, 7f ore , which has to be subtracted. Since /back cannot be di- 
rectly estimated we need a measurable quantity Iq that can be 
used as a proxy for /back- Then the optical depth is given by 



In 



h - h 



(3) 



\I\RDC - lie 

Following lPeretto et al.l d2008l) we assume that 7f ore = I z \, where 
I z \ is the zodiacal light, which is systematically calculated for 
every Spitzer observation and available in the image header. 
For the quantity Iq we used the average emission level from a 
patch of MIR emission in the close vicinity of the actual cloud. 
These emission patches, typically around 1 square arcminute in 
size, were chosen manually in order to capture the characteristic 
emission level for the background approximation of the individ- 
ual clouds and to exclude strong compact emission sources. The 
mean and standard deviation of the emission levels within these 
defined regions have been computed in order to be used as Iq 
in Equ. [3] The standard deviation obtained here is propagated 
through the following steps and is used to give formal errors for 
the derived masses and column densities, as listed in Table [3] 

After optical depth determination, this quantity is converted 
into column densities and finally to masses by using the follow- 
ing equations: 



N H , = 1.086-, 
M = m H A Nh, , 



(4) 
(5) 



where Nn 2 is the H2 molecule column density, me, is the mass of 
one hydrogen molecule, and A is the area per pixel, cr is the ex- 
tinction cross section per hydroge n molecule. According to the 
adopted dust extinction model bv lWeingartner & Draind 12001, 
see below), cr = 4.62 x 10~ 23 cm 2 for the Spitzer/IRAC band 4 
central wavelength of 7.87 jum. 

The derived masses critically depend on the used extinction- 
column density calibration, the distance to the targets, and the 
method to select the relevant extinction regions. As dust m odel 
we used the parametrisation of Weinga rtner & Draind d200 ll) ac- 
cording to their model B with Ry = 5.5. This particular model 
has bee n shown to be relevant for massive star-forming regions, 
e.g., by Indeb etouw et al. ( 2005). It d iffers from the common 



dust models (e.g.. lDraine & Lee 



1984) in that it predicts higher 



extinction cross sections especially in the 4-8 micron wave- 
length region, a relevant point for the Spitzer extinction maps. 
We note in passing that such elevated cross sections are also pre- 
dicted in connection with ice-coated dust grains and especially if 
dust c oagulation processes are involved dOssenkopf & He nning 
1994). The chosen dust model finally relates the optical depth 
(and, equivalently, the extinction magnitude at the used wave- 
length) to the equivalent column density of hydrogen molecules. 
The quadratic dependence of the masses on the distance is in the 
extinction map case the same as for masses derived from the 1 .2 
mm emission data. 



3.3. Comparison of techniques 

Using the data from two different regions of the spectrum and 
two different techniques for estimating IRDC parameters, en- 
ables us to compare the results and to analyse advantages and 
disadvantages of both methods (see Table [l}. The mid-IR data 
have an effective resolution of around 3", which is much higher 
than 24" for our millimeter data. The estimated total masses and 
peak column densities have no dependence on the temperature. 
On the other hand, when using 8 //m data to calculate the opti- 
cal depth and then the column density in the cloud, we have to 
estimate the expected flux from the background indirectly. We 
assume the background flux as the difference between an aver- 
age intensity around the cloud and the foreground intensity, but 
in reality matter hidden by the infrared dark cloud can be inho- 
mogeneous and has either lower or higher intensity. However, 
the logarithm in Eq. © to a certain degree mitigates uncertain- 
ties in the estimation of the intensities necessary for evaluating 
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A R.A. [arc minutes] A RA l ate minutes] 

CENTRE: R.A. 1S h 17 m 26f71 DEC -1 7°u5'06:2 (J2000) CENTRE: R.A. 13 n 37 m 00?37 DEC -62°44 , 04:5 (J2000) 



IRDC 309.134-0.149 2.99 




21 0-1-2 °- 00 21 0-1-2 000 



A R.A. [arc minutes] A RA [ar0 minut8s] 

CENTRE: R.A. 13 B 56™05;59 DEC -62°13'44:6 (J2000) CENTRE: R.A. 14 B ir42'07 DEC -61°29'57:'3 (J2000) 

Fig. 3. Column density maps derived from 8 ^m extinction overlaid with 1.2 mm continuum emission as contours. The scaling is 
indicated in the bar to the right of each image. The contours are 60, 108, 156, 240, 360, 480 mJy beam -1 in all cases. 
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Fig. 4. Column density maps derived from 8 fim extinction overlaid with 1.2 mm continuum emission as contours. The scaling is 
indicated in the bar to the right of each image. The contours are 60, 108, 156, 240, 360, 480 mJy beam -1 in all cases except for 
IRDC 316.45-0.65 which is 84, 120, 156, 240, 360, 480 mJy beam-'. 



8 T. Vasyunina et al.: Physical Properties of Southern Infrared Dark Clouds 

Table 2. Comparison of the different observational techniques for IRDC 18223-3. 



Data 


Resolution 


Column density 


Optical depth 


Contrast 




(arcsec) 


(10 22 cm" 2 ) 


at 8 yum 


(Ibg/IlRDc) 


8 yum GLIMPSE 


3 


2.3 


0.95 


2.5 


1.2 mm IRAM" 


11 


5.9 


2.4 


11 


3.2 mm PdBI b 


5'.' 8 x TA 


45 


19 


1.8e+08 


1.3mmSMA c 


1'.'3 x Y'.A 


93 


40 


1.8e+17 



a [ Beuther et ail d2002h 
b beuther et all (12005) 

c Fallscheer, Beuther, et al. 2009 in preparation 



Eq. (H. 

Simo n et al. (2006a) have used an automated approach where 
they strongly smooth the original (MSX) mid-infrared images 
and basically use these smoothed data as approximation for the 
quantity Iq (cf. Eq. [3]). We tested this approach for one of our 
clouds but found that unless very large smoothing kernels are 
used, the herewith derived intensity contrasts are systematically 
smaller than for the method we employed (Sect. 13.2b . For very 
large smoothing kernels (> 10') it is very difficult to control 
the process. Continuum emission from strong extended emission 
sources in the surroundings might be folded into the area of the 
IRDCs, depending on the individual circumstances. Therefore, 
we refrain from using this smoothing approach in order to have 
better control for the estimation of the quantity Iq. 
Uncertainties in the foreground level estimation also affect the 
resulting masses and column densities. At the moment we take 
into account only the zodiacal light contribution, which typi- 
cally accounts for 15%-20% of the intensity toward the IRDCs. 
An extreme (hypothetical) case may be, that all remaining flux 
received from the IRDC, down to the noise level 7 no j se of the 
GLIMPSE images, is produced by another foreground emission 
contribution (e.g., the PDR surface of the cloud, glowing in PAH 
emission). After removing all this emission down to / no i se , the 
resulting intensity contrast would of course be clearly higher. 
However, considering typical levels for / n0 i se and /Back, our re- 
sults for the column densities would rise just by a factor of 
around 3. This is again due to the alleviating effect of the log- 
arithm in Eq. acting on the intensity contrast. 

While for almost all our clouds peak column densities ex- 
tracted from the GLIMPSE data are slightly higher than the 
values derived from the millimeter data, this difference does 
not correspond to the factor 8 in resolution. Therefore, we 
made a comparison between the peak column densities, ex- 
tracted from different observational data for the p reviously stud- 
ied IRDC 18223-3 dBeufher et alj [20021 120051) . which is lo- 
cated 3. 7 kpc from us and h as a NH3 rotation temperature of 
~ 33 K (Sridh aran et al.l2 005). Fo r this particular clou d, 1.2 mm 
IRAM single-dish observations (Beut heret alj|200"2l) . 8 micron 
Spitzer/GLIMPSE data and interferometric data at 3.2 and 1.3 
mm data obtained with the PdBI and SMA (iBeufher et alj [2005. 
Fallscheer et al. 2009 in prep.) are available, with the spatial res- 
olution ranging from 1 1 arcsec down to 1 .4 arcsec (Tabled. For 
the GLIMPSE data we obtained a column density distribution 
map according to the algorithm described in Sect. 13.21 Together 
with the 1 .2 mm IRAM data as contours, this result is presented 
in Fig. [5] The peak column densities for all the millimeter data 
were calculated using Eq. ©, a dopting the different b eam sizes. 
Always a temperature of 33 K ([Sridharan et al. 2005) has been 
used. The peak flux density for the 3.2 mm PdBI data we took 



IRDC 18223-3 




-4 li 1 I 1 1 1 i I i 1 i 1 I 1 i 1 i I ■ 1 



A R.A. [arc minutes] 
REFERENCE POS. R.A. 18 h 25 m 08?60 DEC -1 2°45'23'. , 5 (J2000) 

Fig. 5. Column density map derived from an 8//m 
Spitzer/GLIMPSE image overlaid with 1 .2 mm IRAM contin- 
uum emission as contours. The image scaling is indicated in the 
sidebar. The contours are 38,76,114,190,266 mJy beam -1 . The 
(0,0) position corresponds to the center of the IRDC 18223-3. 



from the corresponding paper (Beuth er" et alj 12005). while we 
(re-)assessed the peak flux densities for the IRAA/Q and SMA 
data on the related FITS files, kindly provided by H. Beuther and 
C. Fallscheer. For the dust opacity per gram of dust, as for our 
SIMBA millimeter data, we have always used the same opac- 
ity mode l, appropriate for coagulated dust particles with thin ice 
mantles (Ossenkop f & Henningll 19941 opacity column 5 in their 
Table 1). As we can see in Table [2] the peak column density 



2 The 1.2 mm peak flux de nsity we find in the I RAM 30-m data is 
clearly higher than reported in Beuther et al] d2002h . In that paper, the 
millimeter peaks had been fitted with Gaussians which occasionally un- 
derestimated the true peak flux densities. 
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derived from the high-resolution millimeter interferometry data 
is a factor of 40 higher than the one extracted with the mid-IR 
technique. A factor of a few between the mid-infrared and the 
millimeter single-dish results for the peak column densities can 
be accounted for by using other dust opacities/extinction cross 
sections or more extreme assumptions on the MIR foreground 
contribution (see above). However, the large difference between 
the mid-infrared and the millimeter interferometry results indi- 
cates a principle limitation of the extinction method to distin- 
guish high column density peaks. The realistically attainable in- 
tensity contrast at high optical depths sets this limit. In Table [2] 
we list the optical depths and image contrasts at 8 micron that 
would correspond to the column densities derived from the mil- 
limeter data. Obviously, realistic 8 micron images cannot pro- 
vide such humongously high dynamic ranges necessary to derive 
column densities similar to the interferometry results. 

4. Results 

4.1. Morphology of the IRDCs 

The morphology of the clouds in our sample ranged from com- 
pact structures (IRDC 312.33-0.07 PI and P2) to filaments 
(IRDC 309.13-0.15). They have sizes from 1' (IRDC 013.84- 
0.49 P2 and P3) to 4' (IRDC 317.71+0.11), which roughly cor- 
responds to 1-3.5 pc at the distance of these clouds. There 
is, in general, a good agreement between the morphologies of 
the 1.2 mm emission and the 8 /urn extinction structures (see 
Fig. [2H]). As a rule, dense areas of extinction material (IRDCs 
013.84-0.49, 313.72-0.29 etc.) coincide with relatively bright 
compact sources at 1.2 mm. For some IRDC complexes, some 
emission peaks at the millimeter wavelengths are associated 
with mid-IR emission sources (IRDCs 309.13-0.15 P2, 309.34- 
0.64 P2 and 317.71+0.11) which indicates later evolutionary 
stages than those objects corresponding to 8 micron extinction. 
Among all clouds there is one particular case - IRDC 310.39- 
0.28, where the millimeter emission still peaks at the extinc- 
tion maximum despite the existence of a bright MIR emission 
source nearbyQ The IRDCs 013.84-0.49, 312.33-0.07, 318.13- 
0.34, 320.23+0.32, and 321.71+0.07 present several separated 
sources at mm wavelengths, which coincide with dense 8 fim 
features. IRDC 308.12-0.33 has an elongated shape and three 
sub-structures can be recognized, one of them corresponds to the 
extinction maximum at mid-IR wavelengths. IRDC 309.13-0.15 
shows two distinct 1 .2 mm emission sources: the compact object 
(PI) coincides with the MIR extinction region, another one (P2) 
has an emission peak corresponding to the 8 fim emission and an 
elongated "tail". Extended millimeter emission, 2' x 4' in size, 
with at least four sub-structures is present in IRDC 313.72-0.29. 
For IRDC 316.45-0.66, we can recognise only one weak 1.2 mm 
source and diffuse extended MIR extinction structures. 

Positions of known IRAS sources are marked with plus- 
signs on the Figs.[T]-[2] where the IRAC 8 fim data are displayed 
as inverse grey-scale images, and contours present the 1 .2 mm 
data. As a rule IRAS sources do not correspond to the extinc- 
tion regions at 8 fim. On the contrary, they agree with the loca- 
tions where MIR emission coincides with millimeter emission 
peaks or just with the very bright MIR emission sources. In 
the case of IRDC 317.71+0.11, the kinematic distance to one 
of the IRAS sources, located in the center of the 8 fim emis- 
sion, can be derived using the CS (2-1) line velocities reported 

3 Based on the very similar LSR velocities for the MIR source and the 
extinction region, derived from our recent molecular line observations, 
both objects are probably associated. 



by iBronfman et alj ( 1 19961) toward this IRAS source. The v L sr 
velocities of this CS measurement and our HCO + data for the 
neighbouring IRDC are less than 0.5 km/s different. The dis- 
tance to the cloud and the IRAS source is then around 2.9 kpc. 
Hence, assuming that the dark cloud is related with this IRAS 
point source, the infrared luminosity of the compound (IR source 
+ IR DC) is 10 4 L , usi ng the IRAS approximation formula 
from lHenning et al] dl990l) . This is the typical order of magni- 
tude for pre-main-sequence OB stars. 

4.2. Masses and column densities 

Toward all 12 clouds millimeter continuum emission was de- 
tected and column density maps were extracted from 8 fim im- 
ages according to the algorithm described in Sect. 13.21 Figures 
[3]|4]present these column density maps for every region, super- 
imposed with the corresponding 1 .2 mm contours. In Table [3] 
the properties of the IRDCs are compiled: name, right ascen- 
sion, declination, distance, peak flux density and integrated flux 
density, mass and peak column density of the 1 .2 mm sources, 
mass and peak column density of the extinction matter. For every 
cloud in the table, the first line corresponds to the total millime- 
ter flux density as well as total masses of the emission and ex- 
tinction matter. The following lines present data for the separate 
millimeter sub-clumps which are labeled with "P". 

Where appropriate, we distinguished separate sub-clumps at 
1.2 mm. In the Fig. [30 and in Table [3] these separate emission 
sources are named with the designation "P" (e.g. PI, P2). For all 
of them we measured peak flux density and integrated flux den- 
sity, derived masses and column densities according to Equ. Q] 
-12 In the case of IRDCs 309.13-0.15 (P2), 309.34-0.64 (P2) 
and 317.71+0.11 (PI) parts which coincide with strong 8 fim 
emission features were not taken into account for measuring 1 .2 
mm flux densities and hence for estimating masses and column 
densities. The typical range of masses of the separate millimeter 
sources with T = 20 K is 50-1000 M , and the column densities 
range between 0.9 and 4.6 xlO 22 cirT 2 . 

For calculating total masses of the extinction material in in- 
dividual clumps we took regions above 3 sigma into account, 
where sigma is in this case the standard deviation in the full 
extinction maps. The total masses of the IRDCs for extinction 
matter were found to range from 300 to 1700 M and the de- 
rived peak column densities correspond to values from 2. 1 to 5.4 
xlO 22 cirT 2 . 



Table 3. Properties of the IRDCs. For every cloud the first line corresponds to the total mass, the following to the masses of separate millimeter sources. 



Name 


R.A. 


Decl. 


D 


Peak Flux density 


Integrated Flux density 


Mass 1.2mm 


N(1.2 mm) 1 


N (1.2 mm) 2 


Mass (8 fim) 


N (8 ^m) 




(J2000.0) 


(J2000.0) 


(kpc) 


(mJy) 


(Jy) 


(M Q ) 


(10 22 cm- 2 ) 


(10 22 cm" 2 ) 


(M ) 


(10 22 cm" 2 ) 


IRDC 308.12-0.33 


i i n c\ i i 

13 37 01.2 


/'I A A A f\ 

-62 44 40 


4.32 


225 


1.88 


580 


1.7 


60.3 


c -1/^+40 
520 -30 


o o+0 16 
2 - 8 -0.16 


IRDC 309.13-0.15 
PI 


13 45 17.1 


-62 21 57 


3.92 


330 3 
162 


4.39 3 
1.37 


1150 3 
360 


2.6 3 
1.3 


41.9 


1750+ 370 


3 +0 13 

J - -0.13 


P2 








124 


0.96 


250 


0.9 








TT"» 1 — \ /"*■ 1 f\f\ 1/1 f\ /~ A 

IRDC 309.34-0.64 
PI 


ii a o i r\ o 

13 48 39.8 


-62 47 26 


3.46 


504 
216 


3.69 
0.94 


750 
190 


3.9 
1.7 


48.5 


750!™ 


o i +0.24 
-'• 1 -0.24 


P2 








155 


0.81 


170 


1.3 








IRDC 310.39-0.28 


13 56 01.7 


-62 14 27 


4.93 


594 


2.48 


1029 


4.6 


186.1 


1320+*° 


2 6 +016 


IRDC 31 ? 33-0 07 

llvl-/ V Jli.JJ \J .\J 1 


14 11 56 8 


-61 29 25 


4.05 


288 


0.77 


210 


2.3 


76.6 


9QA+I00 


7+0.33 
Z -'-0.16 


PI 








288 


0.46 


130 


2.3 




P2 








92 


0.31 


90 


0.7 








IRDC 313.72-0.29 


14 23 05.4 


-61 14 48 


3.33 


172 


1.98 


370 


1.3 


35.7 


700!'° 


4 1+0.16 
^ -0.16 


IRDC 316.45-0.66 
PI 


1/1 A A c r\ A 

14 44 50.4 


-60 30 54 


3.01 


159 
159 


0.96 
0.58 


150 
90 


1.3 
1.3 


32.3 


a in+80 


o n+0 16 
2 - 9 16 


IRDC 317.71+0.11 
PI 


14 51 07.5 


-59 16 11 


2.9 


954 3 
250 


8.00 3 
1.98 


1150 3 
280 


7.5 3 
2.0 


47.9 


1320!|» 


3 6 +015 


IRDC 318.13-0.34 
PI 


14 55 58.4 


-59 28 31 


2.96 


142 
142 


2.48 
1.02 


370 
150 


1.1 
1.1 


26.9 


680_J«> 


9 i +0.16 
-0.16 


IRDC 320.23+0.32 
PI 


15 07 56.7 


-57 54 27 


1.97 


186 
142 


1.68 
1.00 


110 

70 


1.5 
1.1 


24.6 


600! 32 » 


t 9+0.16 
J ' Z -0.I6 


P2 








186 


0.68 


50 


1.5 








IRDC 321.71+0.07 
PI 


15 18 26.7 


-57 21 56 


2.14 


408 
276 


3.12 
1.69 


240 
130 


3.2 
2.2 


56.9 


460_$ 


r 4+ 0.16 

-0 23 


r>""> 
Fz 








408 


1 A 1 

1.43 


1 10 


3.2 








IRDC 013.84-0.49 
PI 


18 17 21.2 


-17 09 23 


2.66 


438 
214 


9.38 
2.97 


1130 
360 


3.5 
1.6 


77.0 


ii5o_;» 


t q+0.13 
J - J -0.13 


P2 








163 


1.12 


130 


1.3 








P3 








192 


1.45 


170 


1.5 








P4 








438 


3.54 


430 


3.5 









Peak column density per beam 

2 Extrapolated column density, obtained by applying the correction factor explained in Sect. 14.3.21 

3 This mm peak corresponds to a mid-infrared bright source near to the IRDC. 
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4.3. Comparison with results for other cores 

In this chapter we will compare our results with previously ob- 
tained characteristics for low- and high-mass pre-stellar cores. 



4.3.1. Comparison with low-mass cores 

During recent years, extinction mapping has mainly been ap- 
plied along the lines of the near - infrared colour e xcess method 
(e.g., Lombardi & Alvesl l200lt lLombardil I2008I) or classical 
star count techniques in the visible o r near-infrared (e.g., 
Dob ashi et al.ll2005tlFroebrich et al.l2005l) . Such approaches are 
most powerful towards medium-extinction regions. In contrast, 
the extinction map method used in the current paper exploits 
the extinction of mid-infrared extended emission and hence does 
not rely on the identification of stellar background sources. This 
method can peak into cores having column densities of up to 
10 23 cm It has been used previously for the a nalysis of low- 
mass starless cores based on ISOCAM data (Bac mann et alj 
2000) which gives us the opportunity to compare our results. To 
make a fair compariso n, two effects have to be considered, (a) 
iBacmann et al.l d2000t) used the standard iDraine & Leel dl984h 
extinction cross sections. These are more tha n a factor of 3 
smaller than the Weing artner & Draind d200ll) values we use. 
T herefore, we recomput ed the peak column densities reported 
Bacmann et alj ( 2000) for the low-mass cores by adopting the 



m 



Weingartner & Draine (2001) dust extinction model, (b) For the 
low-mass cores, which typically reside at distances of less than 
300 pc, the linear spatial resolution is much better than in the 
case of IRDCs with their distances of 2 - 5 kpc. Therefore, the 
actual column density peaks are much better resolved in the low- 
mass case. To assess the effect of poor spatial resolution on our 
derived peak colum n densities, we used a sy nthetic column den- 
sity map derived in Steinack er et al. l (|2005) for one of the low- 
mass cores of IBacmann et alj d2000l) . namely Rho Oph D (d — 
160 pc). We convolved this map with kernels appropriate to em- 
ulate the much coarser linear resolution toward our IRDC targets 
and computed the ratio of the unsmoothed to the smoothed peak 
column density value. These factors (typically in the range from 
1.5 to 3.5) have been multiplied to the peak column densities we 
have originally derived from the IRDC extinction maps. 
A comparison of the column densities for low-mass cores and 
our IRDCs, after taking into account the just mentioned con- 
siderations (a) and (b), is shown in Fig. [6] Here we can see a 
clear trend for high-mass cores to attain higher column densities 
than the low-mass objects. This qualitative difference shows that 
IRDCs are not just far away Taurus-like clouds, but a distinct 
type of clouds with the potential to form a distinct type of stars 
(see below and Sect. 14.41 ). 

It might be of concern that in step (b) we applied a cor- 
rection factor which increases the column densities with in- 
creasing kinematic distance. To statistically fortify our state- 
ment that Infrared Dark Clouds present "stochastically larger" 
values for column densities than low-mass pre-stellar cores 
we used the (Wilcoxo n)-Mann- Whitney U one-tailed test (e.g., 
IWall & Je nkins 20031). It is a non-parametric test for assessing 
whether two samples of observations come from the same distri- 
bution or not. It works well on a small number of observations 
in one sample. The null hypothesis is that the two samples are 
drawn from a single populati on. For the test we use d low-mass 
col umn density values fro m Bacmann et al.l d2000t) adopted to 
the|Weingartner & Drainel d2001l) dust extinction model and col- 
umn densities for IRDCs without applying the linear spatial res- 
olution correction (b) mentioned above. 



For both populations we computed the Rank Sum within the 
nonparametric Mann- Whitney statistic, usually called U (108 
vs 0). The distribution of the U-statistic under the null hy- 
pothesis is known and can be found in special tables (e.g., 
Siegel & Castellan! 1 1988b . Then, we estimated the probability 
that the values for low-mass objects and IRDCs came from the 
same distribution. The obtained probability lower than 0.005 re- 
jects the null hypothesis and shows that our two samples come 
from different distributions. 



4.3.2. Comparison with high-mass cores 

The high-mass clumps we want to compare have all been ob- 
served at 1.2 mm, comprising the foll owing sample s : high- 
mass starless core candidates (HMSCs) (Sridharan et al. 2005; 
Beut her et alj20 02). Infrared Dark Clouds from Rath borne et al.l 
(2006), and results presented in this paper for our 1.2 mm 
data. To (re-)calculate peak column densities for the first two 
sets of data we used Eq. d2]i as we ll. In the case of HM SCs, 
peak flux densities were taken from Beut her et alj (120021) . dis- 
tance and tempe rature estim ates for all these object s come from 
Srid haran et all d2005l) . In iRafhborne etail d2006l) . for all 38 
IRDCs peak flux densities and distances are presented. For the 
calculation, we assume again 20 K for the temperature, as dust 
opacity x v and gas-to-dust mass ratio we adopt the values 1.0 
cm 2 g -1 and 100, respectively. Hence, we use here the same pa- 
rameters as for the analysis of our 1 .2 mm data, except for the 
different beam size, equal to (ll") 2 in steradians, adapted to the 
IRAM 30-m telescope. 

Still, the measured data are affected by the convolution with 
the observational beam, which results in different linear smooth- 
ing scales for objects at different distances. In order to eliminate 
this smoothing effect, we attempted to extrapolate the measured 
peak column densities per beam back to the true values, assum- 
ing an (analytic) column density profile. But what is an appro- 
priate choice for such a power law? A first idea is to use the data 
at different beam sizes we have for IRDC 18223-3 (Sect. [3~3l . 
However, it turns out that these data sets give no unique trend 
for such a power law. The comparison between the IRAM 30-m 
data and the SMA data speak for a relatively steep power law of 



N 



and the comparison with the PdBI data would even 



indicate a much steeper power law index. The fact that these two 
data sets do not result in a similar power law may be related to 
the different degree these observations are able to recover ex- 
tended emission. Furthermore, at the high spatial resolution of 
the SMA observations, we begin to see a dense rotation struc- 
ture (Fallscheer et al. 2009, in prep.) that is distinct from the 
rest of the surrounding clump. However, we also compared the 
column density values for a few other IRDCs, observed both 
at single-di sh and interferometric resolution (Rath borne et al.l 
20081) . Again, no clear trend for a certain power law range 



at sing 
120061 5 



is obvious. According to John stone et al.l d2003l) . for the filament 
structure of the famous IRDC Gl 1.1 1-0.02 the column density 
profile should be very steep and may reach ~ r~ 3 . Such a sig- 
nificant difference to commonly observed values might be an 
additional feature characterizing massive star forming clumps, 
espec ially if filamentary str uctures are involved. On the other 
hand, Bacman n et alj (|2000) showed that N ~ r 1 can be a rea- 
sonable choice for lower-mass starless cores , and also the single 
dish mapping of young massive clumps by Beut her et alj {2002) 
resulted in less steep power laws for the column density quite 
close to the low-mass core results. Hence, by assuming N ~ r _1 
for the extrapolation to the true column densities will provide us 
with a kind of robust lower limit. 
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Fig. 6. Histogram showing the principle distribu tion of column 
densiti es for a collection of low-mass cores from Bacmann et al. 
(2000) (after taking into account consideration (a) from 
Sect. I4.3.U and for the IRDCs presented here (after taking into 
account consideration (b) from Sect. 14. 3. Ti l. The Mann-Whitney- 
U one-tailed test confirms the clear separation of the two distri- 
butions even when the different spatial resolutions (considera- 
tion (b)) are not taken into account (Sect. |4~3TTT >. 



To recalculate all column density values, an artificial column 
density distribution map was created (by just assuming an ana- 
lytical N ~ r 1 profile) with the peak in the centre of the array 
normalised to 1 .0 and the assumption that 1 pixel corresponds to 
2000 AU (a typical size scale for fragmentation). Then, for ev- 
ery source from the three high-mass samples mentioned above, 
this column density distribution map has been convolved with 
a smoothing Gaussian kernel emulating the effect of observing 
at mm wavelengths with single-dish telescopes (11-24" beam). 
The kernel size (FWHM) is computed by taking the ratio of the 
effective linear resolution (in AU) achieved within the observed 
mm maps and the 2000 AU pixel size from the artificial map. 
After convolving the artificial map with such a Gaussian kernel, 
the peak column densities in these synthetic maps are smaller 
than 1 . The correction factor for every source is then the ratio 
of the original peak column density value (1.0) to the peak col- 
umn density value after beam convolution. These factors are then 
applied to the column density values derived according to the 
Eq. (|2]), which then approximates the true peak column density 
values. These extrapolated values for the new regions reported 
here are also given in column 9 of Table [3] 

Figure [7J summarizes the distribution of the extrapolated 
peak column densities we have finally obtained. The distribu- 
tion of the true peak column densities indicates a similar order 
of magnitude for previously obtained values for the two other 
samples of massive star forming regions and our new SIMBA 
results. Most of the clouds have (extrapolated) column densities 
in the range from 1.0 to 12.0 xlO 23 cirT 2 with few exceptions 
reaching 20.0 - 30.0 xlO 23 cirT 2 . 

4.4. Comparison with theoretical models. 

Infrared Dark Clouds, located several kpc away from us, have 
large masses, volume and column densities, but "high" val- 
ues alone do not guarantee that they are really the progenitors 
of massive stars, and we ne ed additional criteria to est imate 
this possibility. According to Krumholz & McKee (2008) only 
clouds with column densities of at least 1 g cmT 2 , which corre- 




Column density : 



Fig. 7. Distributio n of extrapolat e d col umn densities for 
the HMSCs from ISridharan et alj (|2005), the IRDCs from 
Rath borne et al.l d2006l) . and the IRDCs from the present pa- 
per. Note that we have applied the corrections and extrapolations 
mentioned in Sect. l4.3.2l for all these high-mass cores. 



sponds to 3 xlO 23 cmT 2 in our units, can form massive stars. As 
we can see in Table [3] the direct transformation of observational 
results for IRDCs from our list indicates peak column densities 
one order of magnitude lower than this limit. On the other hand, 
as it was shown in Sect. 14.31 that the directly observed peak col- 
umn densities are still affected by the convolution with the beam 
of the millimeter observations. After taking this into account, the 
derived values for the extrapolated column density rise by a fac- 
tor of 10 or higher (see Fig.fTjl. Thus, we can reach the 3 xlO 23 
cm" 2 threshold in the case of the infrared dark clouds. 

Since we have not only compact sources, but also very 
filamentary structures like IRDC 320.23+0.32 PI, this raises 
the question, w hether it is possible to for m massive stars in 
such structures. Baneriee & Pudritz (2008) showed that fila- 
ments play a dominant role in controlling the physics, accretion 
rate, and angular momentum of the much smaller-scale accre- 
tion disk that forms within such collapsing structures. Large- 
scale filamentary flows sustain a very high accretion rate (M ~ 
10~ 2 M yr~') due to the supersonic gas flow onto the protostel- 
lar disk. These rates are 10 3 times larger than predicted by the 
collapse of the singular isothermal spheres and exeed the accre- 
tion rates necessary to squeeze the radiation field of the newly 
born massive star. Thus, for almost all our clouds we still see the 
potential to form massive stars. 

5. Conclusions 

In this paper, we discuss our progress in understanding proper- 
ties of Infrared Dark Clouds. A set of 12 clouds located in the 
southern hemisphere has been selected from the MSX 8.3 mi- 
cron images. For these clouds 1.2 mm maps were obtained with 
the SIMBA bolometer array at the SEST telescope. GLIMPSE 
mid-infrared images for these regions were retrieved from the 
Spitzer Archive. 

The new sources comprise a variety of IRDC morphologies, 
from compact cores to filamentary shaped ones, from infrared 
quiet examples (no Spitzer 8 fim emission sources) to more ac- 
tive ones. As a rule, our sample shows a good agreement between 
the morphologies of 1 .2 mm emission and 8 fim extinction fea- 
tures. The total masses of the IRDCs were found to range from 
150 to 1150 M (emission data) and from 300 to 1750 M (ex- 
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tinction data). We derived peak column densities between 0.9 
and 4.6 xlO 22 cm" 2 (emission data) and 2.1 and 5.4 xlO 22 cm -2 
(extinction data). 

Since the MIR extinction method has been used previously 
for the analysis of low-mass starless cores we check how our 
findings relate to the published results. To make a fair compari- 
son, we used the same dust model and the same spatial resolution 
in both cases. It is shown, that there is a clear trend for the high 
mass cores to attain higher column densities than the low-mass 
objects. This qualitative difference means that most IRDCs are 
not just far away Taurus-like clouds, but a distinct type of clouds 
with the potential to form a distinct type of stars. A simple sta- 
tistical analysis (the Mann-Whitney-U one-tailed test) confirms 
this statement also when the different spatial resolutions are not 
taken into account and thus the IRDC column densities are un- 
derestimated. 

Using the data from two different regions of the spectrum 
and applying two different techniques for estimating IRDC pa- 
rameters enables us to compare these two methods. On the one 
hand, the extinction technique has some advantages over the mil- 
limeter technique. It is a "cheap" method since GLIMPSE at 8 
yum has covered large parts of the Galactic plane in the 4th and 
1st quadrant. The GLIMPSE data have a better spatial resolution 
than millimeter single-dish data, which reveals the often filigree 
substructures of the clouds. Furthermore, this method does not 
depend on assumptions for the temperature of the IRDCs for 
estimating masses and column densities. However, our compar- 
ison shows, that the extinction method has a principle limitation 
to distinguish very high image contrasts and hence to find high 
column density peaks s> 10 23 cm 2 . Hence with the extinction 
method we can give only a lower limit to the column density 
values, inspite of high r esolution. The limit is around Ay = 75 
mag when applying the Weingart ner & Draind J2001) Ry =5.5 
B extinction law (corresponding to roughly 200 mag w hen fol- 
lowin g the common Ry - 3.1 extinctionlaw reviewed in Mathis 
1990). High-spatial resolution (sub-)millimeter observations are 
hence crucial to assess even higher column density ranges and 
to reveal the actual column density maxima. 

To compare column densities extracted with the emission 
method with previously obtained values for IRDCs and HMPOs, 
we extrapolated them back to the true peak column densities 
by assuming a column density profile ~ r , thus, mitigating 
the spatial resolution differences within the different samples. 
The distribution of the true peak column densities indicates a 
similar order of magnitude for our new SIMBA results and the 
two other samples of massive star forming regions. Moreover, 
the true peak column densities exceed the theoretical limit of 3 
xlO 23 cm" 2 (or 1 g cm" 2 ), recently put forth to distinguish po- 
tentially high-mass star-forming clouds. 

Thus, extracted values for masses and column densities both 
for emission and extinction matter show a clear difference be- 
tween IRDCs and known low-mass pre-stellar cores, and con- 
firm our assumptions that Infrared Dark Clouds can present the 
earliest stages of high-mass star formation. 
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